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Abstract 

In this paper, we study the anomalous flavor changing neutral current Yukawa 
interactions between the top quark, the Higgs boson, and either an up or charm 
quark (tqH, q = u, c). We probe these couplings in e“p —)• —)• i/eHq and the 

channel e“p —>• t'eHb. Both channels are induced by charged current interactions 
through e“p collision at the Large Hadron Electron Collider(LHeC). We study the 
signatures with the Higgs decay modes H —)• 77 , bb and r+T“. Our results show 
that the flavor changing couplings KtqH can be probed down to a value of 0.0162 in 
e“p —)• t'ef t'eHq with H —)• bb at a 14 TeV LHeC with a 150 GeV electron beam 
and 200 fb“^ luminosity. This value of the coupling corresponds to the branching 
ratio Br(t —)■ qH) = 1.34 x 10“^. 
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1 Introduction 


The discovery of the Higgs boson at the Large Hadron Colhder(LHC)[T]|2] is a major 
step towards understanding the electroweak symmetry breaking (EWSB) mechanism and 
marks a new era in particle physics. In order to ultimately establish the nature of the 
Higgs, a precise measurement of the Higgs couplings to fermions and gauge bosons as well 
as the Higgs self-coupling is needed. These precision measurements will be some of the 
most important tasks for experiments at the LHC and the future colliders. According to 
the analyses of the ATLAS and CMS collaborations, the couplings of the Higgs boson have 
been measured with an overall precision of about 15%, which means that there still remains 
some room for the existence new physics. Besides the Higgs boson, the measurement of 
the top quark properties is also important. It is the heaviest known elementary particle 
which makes it an excellent candidate for new physics searches. To probe new physics 
through the Higgs boson and top quark, the top-Higgs Yukawa couplings are of special 
interest since they are sensitive to new flavor dynamics beyond the Standard Model(SM) 
not too far above the electroweak scale. Furthermore, top quark, as the heaviest SM 
fermion, it owns the strongest Yukawa coupling. Among the Higgs couplings to quarks, 
the most promising place to reveal new physics at high energy colliders are processes 
involving top quarks. 

The mass of the top quark is heavier than that of the observed Higgs boson, which 
makes the top quark flavor changing neutral current(FCNC) processes t —)■ qH(q = u, c) 
kinematically accessible. In the SM, processes that are induced by FCNC in top quark 
production or decay are extremely suppressed by the GIashow-IliopouIos-Maiani(G.LM.) 
mechanism[3] according to SM computation, with decay rates of the order of or be¬ 
low. However, new physics scenarios, such as the minimal supersymmetric model (MSSM) 
with/without R-parity Violating[i] [5] p] [7] [5] P] [TU] [TT] . two-Higgs-Doublet Model (2HDM) 
capiiBiiia, Warped Extra Dimensions |16] [T7] . Alternative Left-Right symmetric Mod¬ 
els (ALRM)[T8], Little Higgs with T parity (LHT)|T9], etc, could enhance the FGNG rates 
by several orders of magnitude, thus making them detectable using current experimental 
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data. Therefore, studying the top-Higgs FCNC interactions is important both from a 
theoretical as well as an experimental perspective. 

Up to now, the searches for t qH have been investigated experimentally at the 
LHC which gives the strong limits on the top-Higgs FCNC couplings. Among them, 
the most stringent constraint of Br(t —)■ cH) < 0.56%, Br(t ^ uH) < 0.45% at 95% 
conhdence level (C.L.) was reported by the CMS collaboration from a combination of the 
multilepton channel and the diphoton plus lepton channel [20]. While an upper limit is 
set on the t —?■ cH branching ratio of 0.79% at the 95% conhdence level by the ATLAS 
collaboration [21] [22]. Except for the widely studied t —)■ qH decays, the importance of the 
single top Higgs associated production has been also emphasized in the recent theoretical 
studies especially at the LHC [23] [21] [2S] [26] [27] [2B] [29] [30] [SI] [32] • In our present paper, 
we study the anomalous FCNC Yukawa interactions between the top quark, the Higgs 
boson, and either an up or charm quark (tqH, q = u, c) at the Large Hadron Electron 
Collider(LHeC). The LHeC kinematic range exceeds HERA’s by a factor of about 20, due 
to the combination of a 7 TeV or higher proton beam from the LHC and a new 60 GeV to 
150 GeV electron beam. Its luminosity is projected to be as high as possibly 10^^cm“^s“^, 
with a default design value of 10^^cm“^s“^. This is almost a thousand times higher than 
HERA’s luminosity, which gives the LHeG the potential of a precision measurement Higgs 
production facility and enables a very large variety of new measurements and searches to 
be conducted. Typically we choose two channels to study the anomalous FGNG Yukawa 
interactions at the LHeG. One is the channel e“p Uet i^eHq with q = u, c and the 
other is the channel e“p —?• i/gHb- Both channels are charged current (GG) interaction 
processes induced through e“p collision at the LHeG. 

Our paper is organized as follows: we build the calculation framework in Section 2 
including a brief introduction to the anomalous flavor changing tqH couplings and our 
selected production channels. Section 3 is arranged to present the numerical results as 
well as the signal and background analysis. Typically, the H —> 77 , bb, r+r“ decay modes 
are taken into account. In Section 4 we present bounds on anomalous tqH couplings at 
the future LHeG. Finally we summarize our conclusions in the last section. 
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2 Calculation Framework 


2.1 Flavor Changing tqH Couplings 

Considering the FCNC Yukawa interactions, the SM Lagrangian can be extended simply 
by allowing the following terms, 

= ft^tuntuH + KtcnfcH + h.c., (1) 


where the real parameters KtuH and KtcH denote the flavor changing couplings of Higgs 
to up-type quarks. Now we have mt — mh larger than mc,mu,mb. In addition to the 
usual top decay mode t —)■ W^b, the top quark can also decay into a charm or up quark 
associated with a Higgs boson. Therefore, the total decay width of the top-quark T^ is 

ht = -I-Tt^cH + (2) 


The decay width of the dominant top quark decay mode t —?• W b at the LO and the 
NLO could be found in refjM]. It is given below 
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where ro(t —)■ W“b) = ^^^|VtbP/Sw(3 ~ 2/3^) is the LO decay width and /3w = (1 “ 
m^/m^)5 is the velocity of the W boson in the top quark rest frame. Gp is the fermi 

constant. The t —)■ u(c)H partial decay width is given as[ 
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Here the Higgs boson and the top quark masses are chosen to be mn = 125.7 GeV and 
mt = 173.21 GeV respectively. Similar to the top quark decay, the new interactions 
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affect also the width of the Higgs though the additional decay into an off-shell top that 
subsequently leads to a single W decay of the Higgs, namely H u(c)(t* —)■ Wb) where 
t* denotes off-shell top quark. Therefore we get 

Th = -i- + rH^.u(t*^bw+) + + rH^-c(t*^-bw+) (6) 

where T^'^ is the normal Higgs decay width in SM. While other terms related to the Higgs 
boson three-body decays are numerically estimated following ref [23j . 

The stringent constraints on the anomalous FCNC couplings are set exploiting the ex¬ 
perimental data of the CMS and ATLAS Collaborations|20][2T][22]. Theoretically, many 
other phenomenological studies are performed based on these experimental data. The 
analysis of ref[21] emphasizes the importance of anomalous single top plus Higgs pro¬ 
duction at the LHC deriving the 95% C.L. upper limits Br(t —)■ cH) < 0.15% and 
Br(t —)■ uH) < 0.19%. Ref[25] studies the single top and Higgs associated production 
pp tHj in the presence of top-Higgs FCNC couplings at the LHC, giving the upper 
limits as Br(t —?■ qH) < 0.12%, Br(t —uH) < 0.26% and Br(t —?• cH) < 0.23% with 
an integrated luminosity of 3000 fb“^ at = 14TeV. Ref|26] quotes a 95% C.L. limit 
Sensitivity in the tt —)■ Wb -|- hq —)■ t'z/b -|- ^^(77)q hnal state of Br(t —)■ qH) < 5(2) x 10“^ 
with an integrated luminosity of 300(3000) fb“^ at = 14 TeV. As can be seen the 
upper limits on the flavour changing top quark decays can be signihcantly improved as 
expected at a High Lumi(HL)-LHC. Ref[27] derives model-independent constraints on the 
tcH and tuH couplings that arise from the bounds on hadronic electric dipole moments. 
Refs [28] and [29] study the top quark decay into Higgs boson, a light quark and top Higgs 
associated production including the next-to-leading order QCD effects. Other related 
publications can be found, for example, in refs [30] [SI] [32] [33], etc. 

2.2 The Processes 

Now we turn to study the selected production processes where the effect of the flavor 
changing couplings could be signihcant. 
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2.2.1 e p ^ z/et ^ z^eHq Channel 



Figure 1: Partonic Feynman diagrams for e p —)■ z/gt ^ z/gHq with q=u ,c. Black blobs 
represent the anomalous tqH couplings parameterized by Eq.([T]). 

The hrst channel we will consider is e“p —)■ z/gt —)■ z/gHq production. The parton level 
signal process at the LHeC can be expressed as 

6 (Pl) + b(p2) —)■ Z/g + t —)■ Z^e(P3) + H(p 4) + q(p5) (7) 

with q = u, c and Pi are the four-momentum of initial and hnal particles, respectively. The 
Feynman diagram for the partonic process is depicted in FigJT] The flavor changing vertex 
proportional to the flavor changing coupling zttqu occurs via the single top production with 
its following decay to Higgs plus u or c quark, where this single top quark is induced by 
the collision of b quark from the proton with the W“ boson emission from the electron 
beam. We thus expect the cross sections for these processes to be proportional to czt^qjj 
where c is some related constants. The parent level signal process e“p —)■ z/gHq -|- X, the 
kinematic distributions and integrated cross sections can then be obtained by convoluting 
the parton level process with the parton distribution function (PDF)|36] of quark in the 
proton, 

dcT(e“p z^eHq -|- X) = / dx G^/p (x, /if)d(T(e“b —?• z/gt —)• z^gHq, \/I), (8) 
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where = 2-^xEeEp is the center-of-mass(c.m.) colliding energy and x is the momentum 
fraction of anti-b quark from proton. 

2.2.2 e“p —)■ z/gHb Channel 



Figure 2: Partonic Feynman diagrams for e p —)■ z/eHb. Black blobs represent the 
anomalous tqH couplings parameterized by Eq.([T]). 

The second channel we considered is e“p —)■ i/eHb production. The parton level signal 
process at the LHeC can be expressed as 

e“(pi) + q(p2) ^ i^e(p3) + H(p 4) + b(p5). (9) 

The Feynman diagram for the partonic process is depicted in FigJ2J The FCNC top-Higgs 
Yukawa couplings are deduced from the initial state u(c)-quarks from the proton collision 
with the anti-top quark from the Wtb coupling. Similarly, the parent level signal process 
e“p —)■ i/eHb -|- X is present as 

dcr(e“p —>• i/eHb -|- X) = J ^q/P (x,/^f)d(T(e“q —)■ z/eHb, Vi) (10) 

where q = u, c and VS is again the c.m. colliding energy at the LHeC. 

2.2.3 Charged Current and Neutral Current Production at the LHeC 

The two channels e“p —)■ z/gt —)■ z/eHq and e“p —)■ z/gHb that we have presented are 
charged current (CC) processes where the CC production leads to a top-beauty associated 
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production through W“ boson emission from the initial electron. In addition to CC 
production, the flavor changing Yukawa couplings can also be produced through neutral 
current (NC) productions. In NC it gives rise to pair production of top-antitop quarks 
from a neutral photon/Z boson emission from the initial electron. A comparison of the 
cross sections of these CC and NC production channels including the anomalous FCNC 
top-Higgs Yukawa couplings is presented in TabiH Here the input parameters and the 
very basic kinematical cuts will be presented in our following discussion. 


(A-.A)= (150 IGeV], 14 [TeV]) 

channels 

cr(KtuH = 0.1) [fb] 

e“p —)■ z/gt —)■ z/eHq 

41.64 

e“p —)■ z/eHb 

1.987 

e“p —>■ e“Ht 

0.616 

e“p —)■ e“HqW 

0.901 


Table 1: A comparison of the cross sections of CC and NC production channels including 
the anomalous FCNC top-Higgs Yukawa couplings with KtuU = 0.1. 

From TabiH we see that the largest production is CC e“p —)■ Ugi i^eHq production. 
For KtuH = 0.1 it is more than 10 times larger than the sum of the other channels. Different 
from the CC production which leads to a top-beauty hnal state, the NC production gives 
rise to pair produced top-antitop quarks. The NC productions are small, but still sizeable 
at the LHeC especially when the polarized electron beam is considered. Furthermore, 
in sharp contrast to the LHC, the absence of pile-up and underlying event effects at 
the LHeC, high rates of single anti-top production is expected to providing a better 
insight through these production channels. The rapidity distributions of the Higgs boson 
through different channels are given in FigJS] In our paper, we only consider the CC 
interactions which dominate over all the other production mechanisms. This includes 
e“p —)■ z/et —)■ z^eHq and e“p —)■ z/gHb channels. Looking at Tab{Tl we also hnd that the 
cross section of the former channel is larger than that of the latter one by roughly a factor 








Figure 3: The rapidity distributions of the Higgs boson through different channels in¬ 
cluding e“p —)■ ^'eHq, e“p —)■ r'eHb, e“p —)■ e“Ht and e“p —)■ e“HqW productions. 

of 20. At hrst sight this seems odd, because the transition e“p(b) —)■ z/gt —)■ z^eHq involves 
an (anti)bottom-qnark PDF, while the transition of e“p(q) —?• z/eHb does not. One is 
therefore tempted to think that the cross section of e“p —?• net z/eHq is smaller than that 
of e“p z/gHb. However, this naive assertion is incorrect, because for e“p —?• Ugt —)■ z/gHq 
the internal (anti)top is exchanged in the s-channel, while in the case e“p —)■ z/gHb the 
top appears in a t-channel exchange. The PDF snppression is thus over compensated by 
an on-shell enhancement. 


3 Results 

3.1 Input Parameters 

We take the input parameters as |37] aew(mz)~^lMS ~ 127.9, Gp = 1.166370 x 10“® GeV“^, 
mz = 91.1876 GeV, so we have m^ = 79.82436 GeV and sin^6'w = 1 — (mw/mz)^ = 
0.233699. For the strong conpling constant we take Og = 0.1184. Thronghont this paper, 
we set the qnark masses as mu = md = mg = mg = 0 GeV and mb = 4.18 GeV. The top 
qnark mass is set to mt = 173.21 GeV with its width Ft = 1.3604 GeV when KtuH = 0.1. 
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For the leptons, we keep nie = = 0 GeV, and m,- = 1.77682 GeV. We do not consider 

the contribution from small GKM matrix Vqq/ where q and q' are not the same generation. 
For the mass of the Higgs boson, we take mn = 125.7 GeV with the SM width to be 
F™ = 4.3 MeV. The partonic cross sections are convoluted with GTEQ6L1[3S] parton 
distribution functions (PDF) keeping factorization and renormalization scale /if = /if = 
mf. For the LHeG colliding energy, we consider the future 14 TeV proton at future LHG 
and an energetic new electron beam with the energies of 150 GeV|39][in]. The luminosity 
is taken to be a running parameter. The FGNG couplings are chosen to be /tfuH = 0.1 
and KfcH = 0 for simplicity. This set of parameters will be used as default unless being 
stated otherwise. 

3.2 Kinematic Cuts 

The event reconstruction is still based on a parameterised, generic LHG-style detector. 
The general acceptance cuts in the lab frame for the events are: 

pi^* > 25 GeV, p!^ > 25 GeV, p^ > 25 GeV, p^ > 25 GeV, > 25 GeV, 

< 5, \ri^\ < 2.5, |r/^| < 2.5, \t]^\ < 2.5, 

AR(jj) > 0.4, AR(bb) > 0.4, AR(££) > 0.4, AR(77) > 0.4, AR(7£) > 0.4 
AR(jb) > 0.4, AR(£j) > 0.4, AR(£b) > 0.4, AR(7b) > 0.4, AR(7j) > 0.4 (11) 

where AR = y^A^M^^Ar/^ is the separation in the rapidity-azimuth plane, are the 

transverse momentum of jets (refer as j), leptons and photons while is the missing 
transverse momentum. We stress here that cuts in Eq. fllll) are the very basic ones and 
might be changed later in our following discussion. 
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3.3 Decay Modes and Backgrounds 

3.3.1 e“p —)■ z/et —)■ z^eHq Channel with H —)■ 77 (bb, Decay Modes 

Lets first consider the e“p —)■ z^et z/gHq Channel with H —)■ 77 decay mode. The 
considered signal prodnction can be written as 

e"p -)■ z/gt ^ Z/gHq -)■ z/gqqq ( 12 ) 

with q = n, c. Since in onr calculation we take the anomalous FCNC couplings to be 
«^tuH = 0.1 and KtcH = 0 for simplicity, only q = u contributes. Higgs decays to pairs 
of photons are simulated using MadGraphjlTj where the implementation of the effective 
H 77 interaction is adopted For simplicity one can also multiply the production cross 
sections with the Higgs branching ratio corresponding to the hnal state. As can be seen, in 
this case, the studied topology of our signal gives rise to the jet + + diphoton signature 

characterized by one jet, a missing transverse momentum (^'t) from the undetected neu¬ 
trino and a diphoton signal appearing as a narrow resonance centered around the Higgs 
boson mass. The irreducible background comes from the SM process e“p —)■ z/e 77 q which 
yields the identical final states with the signal. These backgrounds mainly come from the 
production of W boson with double photon production through WW 77 , W 77 couplings 
or through WW —> H —)■ 77 decay associated with jet emission. The others come from jet 
production associated with emission of photons. In order to obtain the anomalous FCNC 
tqH coupling effects, we need to simulate all the signal contributions precisely together 
with these irreducible backgrounds as well as their interference. The total cross section 
for these reactions thus can be split into three contributions 

a = a.Q + + a2KtuH (13) 

where ao is the SM prediction, the term ai linear in zttuH arises from the interference 
between SM and the anomalous amplitudes, whereas the quadratic term a 2 is the self¬ 
interference of the anomalous amplitudes. Potentially reducible backgrounds come from 
various other SM processes that yield different hnal states which are attributed to the 
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jet + + diphoton signature due to a misidentification of one or more of the dual state 

objects. For example, two light jets production with both jets faking a diphoton pair, 
one jet one photon associated production with one jet faking a photon or leptons faking 
photons, etc. The background arising from e“p —)■ z/eZ 7 ee“ 7 j is smaller than 1 percent of 
signal after applying all cuts and taking rejection factors into account. We consider all 
these contributions and take the jet faking a photon rate to be 0 . 001 , the electron faking 
photon rate to be 0.062 113] during data analysis. Although the 77 decay channel has a 
small branching ratio, it has the advantage of good resolution on the 77 resonance and 
is also free from the large QCD backgrounds. Typically, we use a narrow invariant mass 
window — mn] < 5 GeV to further reduce the non-resonant backgrounds as well as 
the jet such that the invariant mass of jyy system to be near mass of the top quark, say, 
mj.y.y belongs to the range [mt — 10, mt -|- 10] GeV. 



[Pb] 

Gross sections for e p —)■ z/gt —)■ z/gHq channeh^tuH = 0.1) 

Decay 

Gut I 

Gut II 

Gut III 

Gut IV 

H —)■ 77 


9.31x10-'^ 

9.26x10-3 

9.21x10-3 

9.04x10-3 

O-B 

2.75 X 10-2 

1.35 X 10-2 

1.02 X 10-3 

5.29 X 10-3 


S/B 

3.39 X 10-3 

6.86 X 10-3 

9.03 X 10-2 

1.71 


CTS 

1.33 X 10-2 

1.33 X 10-2 

1.33 X 10-2 

1.30 X 10-2 

H ^ bb 

2.65 X 10-3 

1.97 X 10-3 




O-B 

6.12 X 10-2 

3.02 X 10-3 


S/B 

5.02 X 10-2 

6.75 X 10-2 

2.17 X 10-3 

4.30 



2.24 X 10-3 

2.23 X 10-3 

2.23 X 10-3 

2.23 X 10-3 

H —)■ T~^T~ 






O-B 

4.93 X 10-2 

1.87 X 10-2 

6.89 X 10-3 

4.20 X 10-^ 


S/B 

4.54 X 10-2 

1.19 X 10-3 

3.24 X 10-3 

5.31 


Table 2: Signal and total Background cross sections for e“p —)■ Ugf channel 

with different decay modes after the application of Gut I-IV. The rejection factors and 
b-tagging effects are taken into account in this table. 
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Figure 4: The signal and total background transverse missing energy(^™^*) distributions, 
transverse momentum distributions and Ai?( 77 ) distributions for e“p —)■ 

after considering Cut I-IV. The anomalous coupling is chosen to be KtuU = 0.1. The 
rejection factors are taken into account. 
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We define some sets of kinematical cuts as bellow: 


• Cut I means the basic cuts present in Eq. ffTT]) : 

• Cut II means the basic cuts plus 25 < < SOOGeV, 25 < pi^* < lOOGeV, 

25 < p:^ < 200GeV; 

• Gut III means Gut II plus requiring the invariant mass of the diphoton pair to be 
in the range [mn — 5, mn + 5 ] GeV; 

• Gut IV means Gut III plus requiring the invariant mass of the diphoton and light 
jet system to lie in the range [mt — 10, mt + 10] GeV. 

In TabJSl we display the signal and the main background cross sections for e“p —?• z/eyyj 
after the application of Gut I-IV. The rejection factors and the b-tagging effects are al¬ 
ready taken into account in this table, where as means the cross section for signal, ctb 
for the background. In FigH] we display the signal’s and the total background’s trans¬ 
verse missing energy(.^™*^) distributions, transverse momentum distributions and 

A/?( 77 ) distributions for e“p —)■ i^e77j in parton level after considering Gut I-IV. The 
anomalous coupling is chosen to be KtqH = 0.1. The rejection factors are taken into ac¬ 
count. We see that the anomalous FGNG tqH couplings can enhance the SM production 
to a level where it can be detectable at future LHeG. By a simple £t we get the hnal cross 
section to be Utotai = 5.10 x 10“® -|- 5.21 x lO^^KtuH + 8.63 x [pbj. 

Now we consider the e“p —)■ z/ef —)■ z/eHq channel with H —)■ bb decay mode. In 
this case, the signal production channel is characterized by a missing energy from the 
undetected neutrino and a bb pair associated with a light jet signal. Still, the bb pair 
signal is appearing as a narrow resonance centered around the Higgs boson mass. The 
main background processes are e“p —)• z/ebbj, bjj, bjj, etc, with light jets faking b jets. 
In our analysis, we assume a b-jet tagging efficiency of Cb = 60% and a corresponding 
mistagging rate of eiight = 1% for light jets (u, d, s quark or gluon) and ec = 10% for 
a c-jet, consistent with typical values assumed by the LHG experiments |44j . For this 
decay mode, we take Gut I, III, IV the same while Gut II to be the basic cuts plus 
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25 < < 400GeV, 25 < < 200GeV, 25 < < 140GeV and AR(bj) < 4. 

For the background production, we also need that the special cut for z/ebjj, Uebcj, etc, 
with the light jets system not belongs to the range [m-w — 10, m-w + 10] GeV. This 
cut will not affect the signal much but it will reduce the background obviously. Finally 
we get the signal and total background to be 13 fb and 3.02 fb, respectively, and we 
get the signal background ratio to be 4.3. The hnal cross section can be written as 
CTtotai = 2.87 X 10"^ + 7.68 x 10~^KtnR + 1-24k2^jj [pbj. 

Finally we consider the decay mode in this production channel. Our results show 
that e“p Uet i^eHq channel with H decay mode can be another good choice. 

With the four lists of cuts, we take Gut I, III, IV the same while Gut II to be the basic cuts 
plus 25 < < 300GeV, 25 < pip < lOOGeV, 25 < p^ < 200GeV and AR(f'j) < 4. The 

total cross section can be parametrised as atotai = 3.96 x 10“^ + 1.30 x 10“^KtuH + 2.13 x 
10“^K^^fj [pbj. The cross sections of the above decay modes are presented in Tabj2]with 
different sets of cuts. We see that the H —)■ 77 decay mode provides the smallest signal 
since the branching ratio of H —> 77 is quite small. By applying the cuts, the background 
can be reduced to the same level. For the H —)■ bb, decay modes, the signal can be 
5 times larger than the backgrounds, thus making the signal over background around 5 
for KtuH equal 0.1. The distributions of the signals and backgrounds are similar to FigJU 
and we therefore do not display them. 

3.3.2 e“p —)■ z/gHb Channel with H —)■ 77 (bb, Decay Modes 

We apply the similar method to e“p —)■ vjlh production Ghannel. However, due to 
the critical large backgrounds, we use much harder cuts instead: For H —)■ 77 decay 
mode. Gut I still means the very basic cuts present in Eq. flTT]) : Gut II means Gut I plus 
< 2.5GeV, Pt > lOOGeV, AR( 7 j) < 4GeV; Gut III means Gut II plus invariance 
mass of diphoton pair belong to [mh — 3, m^ + 3 ] GeV; Gut IV means Gut III plus 
Pt > 150GeV, Pt > 250GeV, AR( 77 ) < 1.5GeV. For H —)■ bb decay mode, we use Gut 

II to be the basic cuts plus < 2.5GeV, AR( 7 j) < 4GeV, and Gut IV to be the Gut 

III plus Pt > 200GeV. For H decay mode, we use Gut II to be the basic cuts 


15 


Decay 

[pb] 

Gross sections for e p —)■ 

Gnt I Gnt H 

Z/gHb channel(fttuH = 0.1) 

Gnt HI Gnt IV 


o-s 

l.bbxlO"*^ 

1.24x10-3 

1.23x10-3 

0.80x10"^ 

H ^ 77 

CTb 

2.70x10"^ 

5.31x10-3 

1.99x10-3 

3.32x10"® 


S/B 

6.15 X 10-3 

2.34 X 10-2 

6.18 X 10-3 

24.1 


o-s 

2.38 X 10-^ 

2.36 X 10"^ 

2.16 X 10"^ 

3.44 X 10-3 

H ^ bb 

CTb 

6.05 X 10-3 

3.14 X 10-3 

1.27 X 10-3 

1.49 X 10"^ 


S/B 

3.93 X 10-2 

7.52 X 10-2 

1.70 X 10-3 

23.1 


c^s 

4.01 X 10-3 

3.99 X 10-3 

4.00 X 10-3 

3.41 X 10-3 

H —)■ r+r" 

ctb 

5.42 X 10"^ 

2.86 X 10"^ 

6.44 X 10-3 

2.23 X 10"^ 


S/B 

7.40 X 10-2 

1.40 X 10-3 

6.21 X 10-3 

15.3 


Table 3: Signal and total Backgronnd cross sections for e“p —)■ z/gHb channel with 
different decay modes after the application of Cnt I-IV. The rejection factors and the 
b-tagging effects are taken into acconnt in this table. 

pins < 2.5GeV, AR( 7 j) < 4GeV, Gnt IV to be the Gnt III pins p^ > 200GeV, 
Pt > 125GeV,AR(f'£) < l.SGeV . When jet fakes b, we replace the cnts for b to jets. 
In Tab|31 we display the signal and total backgronnd cross sections after the application 
of Gnt I-IV. Here in the table the rejection factors and b-jet tagging efficiency are taken 
into acconnt. 

We see that in order to test the anomalons tqH conpling, the best choice of decay 
mode throngh e“p —)■ z/gHb channel is H —)■ bb. Thongh its cross section is mnch smaller 
than that of the one in e“p z/gt —)■ z^eHq channel with associated bb decay mode, 
its signal over backgronnd ratio is not small. However, its cross section is small after 
the critical set of Gnt IV which makes the detection a challenge. By a simple £t we 
get cTtotai = 5.41 X 10"® 7.43 x 10"®KtuH + 8.00 x 10"®Kt^H [pb] for H 77 . atotai = 

1.46 X 10"^ -|- 1.54 X 10"’^fi:tuH + 3.43 x 10""^/?^^^ [pb] for H —bb. atotai = 3.68 x 10"’^ -|- 
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1.86 X 10 ^KtuH + 3.42 X 10 [p^] for H —)■ t^t . 


3.4 Data Analysis and Search Sensitivity 




Figure 5: The contour plots in luminosity-KtuU plane for expected 95% C.L. limits at 14 
TeV LHeC. 


We follow refs [45] [46] exactly to obtain the sensitivity limits. Typically, the limits 
are achieved by assuming the number of observed events equal to the SM background 
prediction, Nobs = (Tb x £ x e, with C for a given integrated luminosity and e the detection 
efficiency, ctb is the cross section of SM background prediction. As can be seen, the SM 
background events can be less or larger than 10 for different values of the luminosity. We 
thus estimate the sensitivity limits on the anomalous tqH coupling through both channels 
by using two different statistical analysis methods depending on the number of observed 
events Nobs- For Nobs < 10, we employ a Poisson distribution method. In this case, the 
upper limits of number of events Nup at the 95% C.L. can be calculated from the formula 

Ppoisson(Nup; k) = 1 - CL. (14) 

Values for limits candidate Nup can be found in Ref.[3Z]- The expected 95% C.L. limits on 
KtqH can then been calculated by the limits of the observed cross sections. The integrated 
luminosity £ will be taken as a running parameter. For Nobs > 10, a chi-square (y^) 
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analysis is performed with the definition 

(15) 

(TbO 

where cxtot is the cross section containing new physics effects and h = is the statistical 
error with N = ctb x £ x e. The parameter sensitivity limits on anomalons tqH conpling 
as a fnnction of the integrated Inminosity can then be obtained. 

In FigJSl we plot the contonrs of expected 95% C.L. limits to KtuH at 14 TeV LHeC 
with 150 GeV electron beam for e“p z/eHq [left panel] and e“p —)■ z/eHb [right 

panel] channels respectively. The solid curve, dotted curve and dashed curve are for 
H —!■ 77 , H bb, H —)■ r+r“ decay modes respectively. From Figl^l we can see that 
the probed /ttuH limits from e“p channel is much smaller than that from 

e“p —)■ i/gHb channel. Typically, we get 0.0588, 0.0162, 0.0209 for Ktuh by using H —)■ 77 , 
H —)■ bb, H —)■ T~^T~ decay modes respectively, which corresponds to the branching ratio 
Br(t ^ qH) = 0.177%, Br(t ^ qH) = 0.0134%, Br(t ^ qH) = 0.0223% at 14 TeV LHeC 
with 200 fb“^ luminosity for former channel, and 0.177, 0.0701, 0.0776 for the latter, 
corresponding to the branching ratio Br(t —?■ qH) = 1.604%, Br(t —?• qH) = 0.252%, 
Br(t qH) = 0.308%. Thus, we apply higher luminosity for the latter channel, see, 
reaching to 1000 fb“^. Then the research limits change to 0.118, 0.0468 and 0.0518 for KtuH, 
which corresponds to 0.713%, 0.112%, 0.137% for the branching ratio. We can see that 
the LHeC sensitivity to the coupling ^tuH is much improved by using e“p —)■ z/gt —)■ z^gHq 
channel. And for different decay modes, H —)■ bb is best one for both channels. 

In TablH we give the Br(t —?■ qH) for different decay modes for both channels at 
14 TeV LHeC with 10(200) fb”^ luminosity respectively. We see that the limits have 
improved by almost 4 times when the luminosity increases from 10 to 200 fb~^. When 
comparing different decay modes, H —)■ bb is the best decay modes for both channels. 
When we come to different channels, e“p —)■ z/gt —)■ z^gHq is much better than e“p —)■ z/gHb 
channels by almost 10 times. Finally, we use our best limits in H —)■ bb decay modes for 
e“p —)■ z/gt —)■ z/gHq channel, and we get 0.0134% for Br(t —)■ qH) as our result at 14 TeV 
LHeC with 200 fb“^ luminosity. 
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e p —)■ z/gt ^ z^eHq 

e 

p —)■ Z/gHb 


£[fb-i] 

Br(t —>■ qH) 

£[fb-‘ 

] Br(t —)■ qH) 

H — > 77 

10 

0.813% 

10 

7.200% 

200 

0.177% 

200 

1.604% 


10 

0.0425% 

10 

1 . 121 % 

H ^bb 

200 

0.0134% 

200 

0.251% 


10 

0.0899% 

10 

1.377% 

H —^ T~^T~ 

200 

0.0223% 

200 

0.312% 


Table 4: Summary for the expected 95% C.L. limits of Br(t —?• qH) for e“p Uet r'eHq 
and e“p —)■ z/gHb channels with H —)■ 77 , H —)■ bb, and H —)■ decay modes at 14 

TeV LHeC with 10(200) fb“luminosity. 

4 Summary and Conclusion 

In this paper, we have investigated the anomalous flavor changing nentral cnrrent (FCNC) 
Yukawa interactions between the top quark, the Higgs boson, and either an up or charm 
quark (tqH, q = u, c). We choose the channel e“p —)■ Ugi —t r'gHq with q = u, c and the 
channel e“p —)■ z/gHb, where both channels are induced by the charged current interaction 
throngh e“p collision at the Large Hadron Electron Collider (LHeC). We consider the 
H —>■ 77 , bb and decay modes. From the resnlts, we can see that the flavor changing 

couplings KtuH can be probed to be minimal as 0.0162(0.0136) for the 95% C.L. limits 
in the e“p Ugt —?• r^gHq channel with H ^ bb decay mode, which corresponds to the 
branching ratios Br(t —)■ qh) = 1.34(0.947) x 10“^ at 14 TeV LHeC with 200(3000) fb“^ 
luminosity. From CMS and ATLAS Collaborations, we get the most stringent constraint of 
Br(t —)■ cH) < 0.56%, Br(t —)■ uH) < 0.45% at 95% conhdence level (C.L.) [20] • Thus, we 
can see that our results shows a strong (above 30 times) improvement from experiments. 
When comparing with the other phenomenological stndies, we can see that the LHeC 
sensitivity our results for Br(t —?• qH) is smaller than the sensitivity limits of LHC as 


19 









Br(t —)■ qH) < 5(2) x 10“'^ with an integrated luminosity of 300(3000) fb“^ at = 
14TeV[26]. Furthermore, our results are comparable with those of other studies, such as 
refs [23] [25]. For example, ref m obtains the sensitivity bound of about 0.1—0.3% through 
different search channels for an integrated luminosity of lOOfb”^ at the ^/s = 13 TeV LHC 
data. 
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